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Abstract 

We present the microscopic equation of growing interface with quenched 

noise for the Tang and Leschhorn model [L. H. Tang and H. Leschhorn, Phys. 

Rev. A 45, R8309 (1992)]. The evolution equation for the height, the mean 

height, and the roughness are reached in a simple way. An equation for 

the interface activity density (or free sites density) as function of time is 

obtained. The microscopic equation allows us to express these equations into 

two contributions: the diffusion and the substratum contributions. All these 

equations shows the strong interplay between the diffusion and the substratum 

contribution in the dynamics. 
PACS numbers: 47.55. Mh, 68.35.Fx 
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I. INTRODUCTION 



The investigation of rough surfaces and interfaces has attracted much attention, for 
decades, due to its importance in many fields, such as the motion of liquids in porous media, 
growth of bacterial colonies, crystal growth, etc. Much effort has been done in understanding 
the properties in these processes M. When a fluid wet a porous medium, a nonequilibrium 
self-affine rough interface is generated. The interface has been characterized through scaling 
of the interfacial width w = ([hi — (hi)] 2 ) 1 ! 2 with time t and lateral size L. The result is 
the determination of two exponents, f3 and a called dynamical and roughness exponents 
respectively. The interfacial width w ~ L a for t 3> L a ^ and w ~ t for t < L a IP. The 
crossover time between this two regimes is of the order of L a ^ . 

The formation of interfaces is determinated by several factors, it is very difficult to 
theoretically discriminate all of them. An understanding of the dynamical nonlinearities, 
the disorder of the media, and the theoretical model representing experimental results is 
difficult to arrive at due the complex nature of the growth. The disorder affects the motion 
of the interface and leads to its roughness. Two main kinds of disorder have been proposed: 
the "annealed" noise that depends only of time and the "quenched" disorder due to the 
inhomogeneity of the media in which the moving phase is propagating. Some experiments 
such as the growth of bacterial colonies and the motion of liquids in porous media, where 
the disorder is quenched, are well described by the directed percolation depining model. 
This model was proposed simultaneously by Tang and Leschhorn || and Buldyrev et al. |§. 
Braunstein and Buceta @ showed that the power law scaling for the roughness only holds 
at criticality for t <C L. Also, starting from the macroscopic equation for the roughness the 
dynamical exponent has been theoretically calculated. They found j3 = 0.629 for the critical 
value q c = 0.539. 

In this paper, we use the TL model in order to investigate the imbibition of a viscous 
fluid in a porous media driven by capillary forces. We write a microscopic equation (ME), 
starting from the microscopic rules, for the evolution of the fluid height as function of time. 
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The ME allows us tO identify two contributions that dominates the dynamics of the system, 
the "diffusion" and the "substratum" contributions. In this context we study the mean 
height speed (MHS), the interface activity density (IAD), i.e the density of actives sites of 
the interface, and the roughness as function of time. We show that the diffusion contribution 
smooth out the surface for q well below the criticality but enhances the roughness near the 
critical value. To our knowledge, the separation into two contributions for all the quantities 
studied in this paper and the important role of the diffusion contribution to the critical 
power-law behaviour has never been studied before. 

The paper is organized as follows. In section |lj we derive the microscopic equation for 
the evolution of height for the TL model. In section |T| we separate two contributions of the 
MHS: the diffusion and the substratum one. We find a relation between these contributions 



that allows us to write an analytical equation for the IAD. In Section [Ty| the temporal 
derivative of square interface width as function of time is derived from the ME and the two 
contributions are identified. These two contributions allow us to explain the mechanism of 
roughness. Finally, we conclude with a discussion in Section |V[ 



II. THE MICROSCOPIC MODEL 

In the model introduced by Tang and Leschhorn (TL) || the interface growth takes place 
in a square lattice of edge L with periodic boundary conditions. We assign a random pinning 
force g(r) uniformly distributed in the interval [0, 1] to every cell of the square lattice. For 
a given applied pressure p > , we can divide the cells into two groups: those with g(r) < p 
(free or active cells), and those with g(r) > p (blocked or inactive cells). Denoting by q the 
density of inactive cells on the lattice, we have q = 1 — p for < p < 1 and q = for p > 1. 
The interface is specified completely by a set of integer column heights hi (i = 1,...,L). 
At t — all columns are assume to have the same height, equal to zero. During growth, a 
column is selected at random, say column i, and compared its height with those of neighbor 
columns (i — 1) and (i + 1). The growth event is defined as follow. If hi is greater than either 
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hi-i or hi+i by two or more units, the height of the lower of the two columns (i — 1) and (i+1) 
is incremented in one (in case of the two being equal, one of the two is chosen with equal 
probability). In the opposite case, hi < min(7ij_i, h i+ i) + 2, the column % advances by one 
unit provided that the cell to be occupied is an active cell. Otherwise no growth takes place. 
In this model, the time unit is defined as one growth attempt. In numerical simulations at 
each growth attempt the time t is increased by St, where St = 1/L. Thus, after L growth 
attempts the time is increased in one unit. In our simulations we used L — 8192 and a time 
interval much less than the crossover time to the static regime. 

We consider the evolution of the height of the i-th site for the process described above. 
We assume periodic boundary conditions in a one-dimensional lattice of L sites. At the time 
t a site j is chosen at random with probability 1/L. Let us denote by hi(t) the height of the 
i-th. generic site at time t. The set of {hi, i = 1, . . . , L} defines the interface between wet 
and dry cells. The time evolution for the interface in a time step St — 1/L is 

hi(t + St) = hi{t) + -Gi(hi-!, hi, h i+l ) , (1) 

where 

G i = W i+1 + W i - 1 + Fi(tijWi , (2) 

with 

W i±1 = G(hi ±1 -hi- 2){[1 - Q(hi - h i±2 )\ + S huhi±2 /2} , 
Wt = l- Q(hi - min(^_ 1 , h i+1 ) - 2) . 

Here h\ — hi + 1 and 6 (re) is the unit step function defined as 6 (x) = 1 for x > and 
equals to otherwise. F^h'j) equals to 1 if the cell at the height h\ is active (i.e. the growth 
may occur at the next step) or if the cell is inactive. Fj is the interface activity function. 
Gi takes into account all the possible ways the site i can grow. The height in the site i is 
increased by one with probability 

1. 1 it j — i + 1 and hi+i > hi + 2 and hi < h i+2 , 
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1/2 



if 3 



% + 1 and hi + i > hi + 2 and hi = /i«+2, 
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1 



if 3 



i — 1 and > hi + 2 and ^ < /i i _ 2 , 
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1/2 



if 3 



i — 1 and > hi + 2 and /ij = /i 



j-2) 
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1 



if 3 



i and /ij < min(/ij_i, /ij+i) + 2 and Fii^h'^ = 1. 



Otherwise, the height is not increased. The cases (l)-(4) are related to growth due to the 
neighbors of the site i. We shall call these mechanisms, growth by "diffusion". Notice that 
these growths are not related to the disorder of the substratum. The factor 1/2 takes into 
account the tie of first-neighbor heights at the (i ± l)-th site in the cases (2) and (4). The 
case (5) is related to local growth, i.e., if the site i is chosen and the difference of heights 
between the i-th and the lowest of his neighbors is less than two, then the height of the 
chosen site increases by one provided that the cell above the interface is active. We shall 
call this mechanism, growth by "substratum" . 

III. MEAN HEIGHT SPEED AND INTERFACE ACTIVITY DENSITY 

Replacing L = 1/St and taking the limit St — > 0, Eq. ([[]) becomes dhi/dt = Gi. Averaging 
over the lattice we obtain (h = {hi)) 



This equation allow us to identify the of two separate contributions: diffusion (1 — Wj) and 
substratum (F{ Wi) 0. Yang and Hu || defined two kinds of growth events: an event in 
which the growth occurs at the chosen site (type A-defined by us as substratum growth) and 
the event in which the growth occurs at the adjoint site (type 5-our growth by diffusion). 
They counted, in numerical simulation, the events number N^it) of type A and A^^(t) of type 
B, in a time interval L. They did not identify this terms as contributions to the mean height 
speed (MHS). Notice that N A (t) oc (Fj W t ) and N B (t) oc (1 - Wi) (see Figure |T|). We shall 
see in Section |IV| that the separation of those two terms allows us to show how the diffusion 



dh 
~di 



(1 - + (FiWi) . 



(3) 
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enhances the roughness near the critical value. The separation into two contributions for all 
the quantities studied in this paper has never been done before. 

The substratum contribution can be expressed as / — (Fi (1 — Wi)), where / = (Fi) is 
the IAD. We found an amazing numerical result: 

(F i (l-W l ))=p(l-W l ) . (4) 

We could not analytically obtain this result. Notice that F{ and 1 — Wi are not independent, 
and that / ^ p for t > 0, as we shall see bellow. Using the Eq. (f|), the IAD is 

f = p(l-W i } + (F l W t ) . (5) 

Figure |2| shows both sides of this equation as function of time showing that Eq. (|j) holds. 
Notice the similarity between Eq. (|3|) and Eq. (j5|). Figure [I] shows the diffusion and the 
substratum contributions as a function of the time for various values of q. At the initial 
time dh/dt — f — p. In the early time regime the substratum contribution dominates 
the diffusion one, because 1 — Wi is very small. The substratum contribution dominates 
the behavior of / and dh/dt in the early regime. As growth continues, the probability 
that growth will occur by diffusion becomes larger; the diffusion contribution increases and 
the substratum one decreases. This can be explained heuristically: inactive sites generate 
a difference of heights greater than two between any site and his neighbor, enhancing the 
growth by diffusion. As time goes on, long chains of pinned sites are generated, slowing down 
the diffusion contribution and hence the substratum one. For q < q c these contributions, 
which in turn dominate, saturate to equilibrium in the asymptotic regime; while, for q > q c , 
both contributions go to zero because the system becomes pinned. At the critical value both 
contributions gives rise to a power law in the IAD and the MHS. Notice that only at the 
critical value does a power-law scaling holds for the MHS (see Figure |3]), which contradicts 
0. This was shown for the roughness by Braunstein and Buceta 
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IV. ROUGHNESS 



From the Eq. (|T|), the temporal derivative of the square interface width (DSIW) is: 

— 2 {{h, — (hi)) Gi) . (6) 

Replacing Gi from Eq. (0), the DSIW can also be expressed by means of substratum and 
diffusion additive contributions. The diffusion contribution is 

2 [ ((1 - Wi) mmfc-x, h i+1 )) - (1 - Wi) (K) ] , (7) 

and the substratum contribution is 

2[(hiF i W i )-(h i )(F i W i )}, (8) 

where the relation Q(x — x') + Q(x' — x) — 8{x — x') = l has been used to derive the diffusion 
contribution. In Figure ^ we plot both contributions as a function of time for various values 
of q. At short times, the diffusion process is unimportant because Ah is mostly less than 
two. As t increases, the behaviour of this contribution depends on q. Notice, from Eq. (|7J), 
that the diffusion contribution may be either negative or positive. The negative contribution 
tends to smooth out the surface. Figure [| shows that this case dominates for small q. The 
positive diffusion contribution enhances the roughness. This last effect is very important at 
the critical value. At this value, the substratum contribution is practically constant, but the 
diffusion contribution is very strong, enhancing the roughness. This last contribution has 
important duties on the power-law behaviour. We think that it is amazing how the diffusion 
plays a dominant role in roughening the surface. To our knowledge the strong effect on the 
roughness, at the criticallity, of the diffusion contribution has never been proven before. 

Generally speaking, the substratum roughens the interface while the diffusion flattens it 
for small q, but the diffusion also roughens the interface when q increases. The diffusion is 
enhanced by substratum growth. The growth by diffusion may also increase the probability 
of substratum growth. This crossing interaction mechanism makes the growth by diffusion 
dominant near the criticality. 
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V. CONCLUSIONS 



We wrote the ME for the evolution of the height in the TL model. The ME allows 
us to separate the substratum and the diffusion contributions and to explain the great 
interplay between them. We found that both contributions to the MHS are related in 
simple way. We found an amazing numerical result that allows to derive the IAD in a 
simple way. The analytical prove of this numerical result is still open. All the quantities 
studied shows, the strong interplay of the diffusion and the substratum contribution in the 
dynamics. The substratum growth enhances the diffusion; increasing the growth by diffusion 
may increase the probability of substratum growth, and vice versa. This crossing interaction 
mechanism makes both dominant at the criticality. The diffusion contribution of the DSIW 
shows different behaviour depending of q. In the intermediate regime, when q is small, this 
contribution is negative, smoothing out the surface. It is astonishing that as q increases the 
contribution became positive, roughing the surface. Finally, we are sure that other DPD 
growth models would permit separation in two contributions with the same features of the 
TL model. 
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FIGURES 

FIG. 1. ln-ln plot of (Ft Wi) (□) and (1 - Wi) (O) versus t. The parameter g is (A) 0.51 (B) 
0.539 (C) 0.6. 

FIG. 2. ln-ln plot of f/p versus t. The symbols show the right-side of Eq. (4) (in units of 
p) compute as is explained in reference [5]. The full curve shows the left-hand side of the same 
equation (in units of p) where / = (Fj). The parameter q is 0.51 (□), 0.539 (O) and 0.6 (y). The 
critical case shows that the IAD goes as t -T? with rj ~ 0.40. 

FIG. 3. ln-ln plot of p^ 1 dh/dt versus t. The parameter q is 0.51 (A), 0.539 (O) and 0.6 (y). 
All cases shows the same behavior in the early time regime. The subcritical case shows that the 
MHS asymptotically goes to certain constant. The critical case shows that the mean height goes 
as . The supercritical case shows that the mean height is asymptotically constant. 

FIG. 4. DSIW (full curve), and its diffusion (O) and substratum (□) contributions versus hit; 
for q equal to 0.3 (A), 0.539 (B) and 0.6 (C). 
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